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Abstract 

We discuss some applications of the Mittag-Leffler function and related 
probability distributions in the theory of renewal processes and continuous 
time random walks. In particular we show the asymptotic (long time) 
equivalence of a generic power law waiting time to the Mittag-Leffler waiting 
time distribution via rescaling and respeeding the clock of time. By a second 
respeeding (by rescaling the spatial variable) we obtain the diffusion limit 
of the continuous time random walk under power law regimes in time and 
in space. Finally, we exhibit the time-fractional drift process as a diffusion 
limit of the fractional Poisson process and as a subordinator for space-time 
fractional diffusion. 
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1. Introduction 

At begin of the past century Gosta Magnus Mittag-Leffler introduced the 
entire function 

oo n 

r(l + an 

n=0 v 

and investigated its basic properties, see Mittag-Leffler (1903). Although 
this function, named after him, and its generalization 

E *Az)--=J2w7T7 w zeC ' 5R« > , > . (1.2) 

T(p + an) 



n=0 



then was investigated by some authors, e.g. Wiman(1905), and used for the 
solution of the second kind Abel integral equation, see Hille and Tamarkin 
(1930), it did not find the deserved attention of the general community. In 
several important books and collections of formulas like Abramowitz and 
Stegun (1965) and Gradshteyn and Ryzhik (2000) on special functions it 
was ignored; a noteworthy exception is Chapter XVIII in Vol. Ill of Erdelyi 
et al. (Bateman project of 1955). 

In this report a prominent role as waiting time distribution will be played 
by the Mittag-Leffler probability distribution function, see Pillai (1990), 

$tf L (t) = 1 - Epi-t?) , f>0, < /? < 1 . (1.3) 

the probability of a waiting time greater than t, called survival probability, 

9f L = l-^ L (t) = E p (-tP), 

and the corresponding probability density function 

<t>f\t) = f*¥ L (t) = ~E p {-iP) = t^ EpA-t") , * > . (1.4) 

The functions Ep(—t") and 0^ L (t) are completely monotone (representable 
as Laplace transforms of non-negative measures, see Feller (1971), in 
concreto Gorenflo and Mainardi (1997), 

i ml/ \ 1 f°° rl3 sin(/37r) . . , , . 

^ ML (t) = - — yH > exp (-rt )dr. 1.5 

^ W vr J r 2 ? + 2r*» cos (/3tt) + 1 P V ; V ; 
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Particularly important is the power law asymptotics for t — > oo 

F(tf>\ ^ rh ML (t\ r(/3 + l) sin(/37r) _ t 
hp{-tr) ~ r _ ^ , (pp [t) ~ (1-0) 



in contrast to the exponential decay of Ei(i) = exp (— t). These asymptotics 
are an essential reason for the importance of these functions in modelling 
anomalous diffusion processes. 

For later use let us here note the Laplace transforms (for < j3 < 1 and 
Res > 0) 

lT ,ML ft i ML 1 

/+oo 
f(x) exp (inx) dx, with /t real, for the Fourier 
-oo 

/■oo 

transform, g(s) := / (/(t) exp (— st) dt, with s in a suitable right half-plane, 
Jo 

for the Laplace transform. 

Remark: Formula (1.5) exhibits for < (3 < 1 the Mittag-Lemer 
waiting time density as a mixture of infinitely many exponential waiting 
time densities rexp(— rt) with r-dependent weight function behaving like 
r~^~ x sin( / 97r) for large r, like r 13 " 1 sin(/37r) for small r. Again we have power 
law asymptotics. For interpretation consider the fact that the density 
rexp(— rt) whose mean is 1/r decays exponentially fast for large r but 
not so fast for small r. Starting around 1965 these functions attained 
increasing attention among researchers, first in the theory of elasticity 
and relaxation (Caputo and Mainardi (1971), Nonnenmacher (1991)), and 
later in the theory of continuous time random walk (pioneered by Montroll 
and Weiss (1965) who, however did not see the relevance of the Mittag- 
Leffler function) and its limiting relation to fractional diffusion. There are 
instances where researchers found the Laplace transform of solutions to 
certain problems but did not identify it as the transforms of functions of 
Mittag-LefHer type, for example Gnedenko and Kovalenko (1968) in the 
theory of thinning or rarefaction of a renewal process, Balakrishnan (1985) 
in his asymptotic investigation of continuous time random walks. For the 
latter it was Hilfer with Anton (1995) who clarified the relationship between 
continuous time random walk (in the sense of Montroll and Weiss), Mittag- 
LefHer waiting time and fractional derivative in time. As more recent 
monographs with useful information on Mittag-Lefher functions let us cite 
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Samko, Kilbas and Marichev (1993), Podlubny (1999), Miller and Ross 
(1993), Kilbas, Srivastava and Trujillo (2006), Mathai and Haubold (2008). 
See also the comprehensive recent report by Haubold, Mathai and Saxena 
(2009). Due to the growing importance of Mittag-Leffler functions there now 
is also activity in the development of efficient methods for their numerical 
calculation, see. e.g. Gorenflo, Loutschko and Luchko (2002) and Seybold 
and Hilfer (2008). 

In Section 2 of the present paper we will sketch the basic formalism of 
continuous time random walk, then in Section 3 under power law regime 
the well-scaled transition to the diffusion limit yielding the Cauchy problem 
for space-time fractional diffusion. Section 4 is devoted to thinning (or 
rarefaction) of a renewal process under power law regime and the relevant 
scaled transition via rescaling and respeeding, to a renewal process with 
Mittag-Leffler waiting time. Then, in Section 5, the Mittag-Leffler waiting 
time law and its relevance in continuous time random walk and the limiting 
fractional diffusion processes are discussed. Again, the transitions are 
achieved via re-scaling and re-speeding. Finally, in Section 6, we discuss the 
time- fractional drift and its role as a time-changing subordinator (producing 
the operational time from the physical time) in space-time fractional 
diffusion. Conclusions are drawn in Section 7. 

2. Continuous Time Random Walk 

Starting in the Sixties and Seventies of the past century the concept of 
continuous time random walk, CTRW, became popular in physics as a 
rather general (microscopic) model for diffusion processes. Let us just cite 
Montroll and Weiss (1965), Montroll and Scher (1973), and the monograph 
of Weiss (1994). Mathematically, a CTRW is a compound renewal process or 
a renewal process with rewards or a random-walk subordinated to a renewal 
process, and has been treated as such by Cox (1967). It is generated by a 
sequence of independent identically distributed (iid) positive waiting times 
T 1; T 2 , T 3 , . . . , each having the same probability distribution <&(£), t > 0, 
of waiting times T, and a sequence of iid random jumps X\, X2, X3, . . . , 
each having the same probability distribution function W(x), x G IR , 
of jumps X. These distributions are assumed to be independent of each 
other. Allowing generalized functions (that are interpretable as measures) 
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in the sense of Gelfand and Shilov (1964) we have corresponding probability 
densities <fi(t) = $'(£) and w(x) = W'(x) that we will use for ease of 
notation. Setting £ = 0, t n = Ti + T 2 + ■ ■ ■ + T n for n G IN, and x = 0, 
x n = X 1 + X 2 + ■ ■ ■ + X n , x(t) = x n for t n <t < t n+ i we get a (microscopic) 
model of a diffusion process. A wandering particle starts in xo = and 
makes a jump X n at each instant t n . Natural probabilistic reasoning then 
leads us to the integral equation of continuous time random walk for the 
probability density p(x, t) of the particle being in position x at instant t > 0: 

ft r r+oo 

p(x, t) = 8(x) V(t) + / 4>{t - t') / w(x-x')p(x',t')dx' dt' , (2.1) 

JO U-oo 

Here the survival probability 

/oo 
4){t') dt 1 (2.2) 

denotes the probability that at instant t the particle still is sitting in its 
initial position xq = 0. Using, generically, for the Laplace transform of a 
function f(t), t > 0, and the Fourier transform of a function g(x), x G H, 
the notations 



POO 

f(s):= / e~ st f(t)dt, 3?s>3?s , 

/+oo 
e +KI 5(x)di, k G 1R , 
-oo 



(2.3) 



we arrive, via 5(k) = 1 and the convolution theorems, in the transform 
domain at the equation 

p(K, s) =$(«) + w(k) 0(s) p(k, s) , (2.4) 

which, by \I/(s) = (1 — 0(s))/s implies the MontroU- Weiss equation, see 
MontroU and Weiss (1965), Weiss (1994), 

~, 1 - d>(s)) 1 , . 

P (k,s) = . (2.5) 

S 1 — w(/t) 0(s) 

Because of |w(/t)| < 1, |0(s)| < 1 for k ^ 0, s ^ 0, we can expand into a 
geometric series 



f( K ,s) = *(s) [<P(s)w(k)Y , (2.6) 



ra=0 
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and promptly obtain the series representation of the CTRW, see Cox (1967), 
Weiss (1994), 

oo 
ra=0 

Here the functions v n — (^ * <p* n ) and w n = w* are obtained by iterated 
convolutions in time t and in space x, respectively, in particular we have 

v (t) = (#*<$)(<) = = (**0)(*)> w (x) = 5(x) , ^(x) = . 

The representation (2.7) can be found without the detour over (2.5) by 
direct probabilistic treatment. It exhibits the CTRW as a subordination of 
a random walk to a renewal process. 

Note that in the special case <f>(t) = mexp(— rut), m > 0, the 
equation (2.1) describes the compound Poisson process. It reduces after 
some manipulations (best carried out in the transform domain) to the 
Kolmogorov- Feller equation 



d_ 

di 1 



/+oo 
w(x — x')p(x', t) dx' 

and from (2.7) we obtain the series representation 

p(x,0=e-™*f;^«; B (x). 



n=0 



3. Power Laws and Well- Scaled Passage to 
the Diffusion Limit 

In recent decades power laws in physical (and also economical and other) 
processes and situations have become increasingly popular for modelling 
slow (in contrast to fast, mostly exponential) decay at infinity. See Newman 
(2005) for a general introduction to this concept. For our purpose let 
us assume that the the distribution of jumps is symmetric, and that the 
distribution of jumps, likewise that of waiting times, either has finite second 
or first moment, respectively, or decays near infinity like a power with 
exponent —a or —j3, respectively, 0<a<2, < /? < 1. Then we 
can state Lemma 1 and Lemma 2. These lemmata and more general ones 
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(e.g. with slowly varying decorations of the power laws (a) and (b)) can be 
distilled from the Gnedenko theorem on the domains of attraction of stable 
probability laws (see Gnedenko and Kolmogorov (1954)), also from Chapter 
9 of Bingham, Goldie and Teugels (1967). For wide generalizations (to 
several space dimensions and to anisostropy) see Meerschaert and Schemer 
(2001). They can also be modified to cover the special case of smooth 
densities w(x) and <p(t) and to the case of fully discrete random walks, see 
Gorenflo and Abdel-Rehim (2004), Gorenflo and Vivoli (2003). For proofs 
see also Gorenflo and Mainardi (2009). 
Lemma 1 (for the jump distribution): 

Assume W(x) increasing, W(— oo) = 0, W(oo) = 1, and symmetry 
W(— x) + W(x) = 1 for all continuity points x ofW(x), and assume (a) or 
(b): 

/+oo 
x 2 dW(x) < oo, labelled as a = 2 ; 
-oo 

/•oo 

(b) / dW(x') ~ ba- x x- a forrr^oo, < a < 2 , 6 > . 

J X 

Then, with \i = a 2 /2 in case (a), p = bn/[T(a + 1) sin(o;7r/2)] in case (b) 

we have the asymptotics 1 — w(k) ~ p\K>\ a for k — > 0. 

Lemma 2 (for the waiting time distribution): 

Assume $(£) increasing, $(0) = 0), $(oo) = 1, and (A) or (B). 

POO 

(A) p :— / td$(t) < oo, labelled as (3 = 1, 

Jo 

(B) 1 - ~ cp-H- p fort^oo, < p < 1 , c> 0. 

Then, with X — pin case (A), A = cn/[T(f3 + 1) sin(/?7r)] in case (B) 
we have the asymptotics 1 — 0(s) ~ Xs 13 for < s — > 0. 

We will now outline the well-scaled passage to the diffusion limit by 
which, via rescaling space and time in a combined way, we will arrive at the 
Cauchy problem for the space-time fractional diffusion equation. Assuming 
the conditions of the two lemmata fulfilled, we carry out this passage in 
the Fourier-Laplace domain. For rescaling we multiply the jumps and the 
waiting times by positive factors h and r and so obtain a random walk 
x n (h) = (X 1 + X 2 + • — h X n ) h with jump instants t n (h) = (7\ + T 2 H — • + 
T n ) t. We study this rescaled random walk under the intention to send h and 
t towards 0. Physically, we change the units of measurement from 1 to 1/h 
in space, from 1 to 1/r in time, respectively, making intervals of moderate 
size numerically small, and intervals of large size numerically of moderate 
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size, in this way turning from the microscopic to the macroscopic view. 
Noting the densities Wh(x) = w(x/h)/h and r (t/r)/r of the reduced jumps 
and waiting times, we get the corresponding transforms Wh{k) = {nh), 
</> T (s) = 0(ts), and in analogy to the Montroll- Weiss equation (2.5) the 
result 

g h|T («, 8 ) = i-M - 1 „ = j-M - 1 „ . (3.1) 
s 1 — Wh(n) 0t(s) s 1 — w(hn) (f)(rs) 

Fixing now /c and s both as 7^ 0, replacing /c by hn and s by rs in Lemma 
1 and Lemma 2, sending h and r to zero, we obtain by a trivial calculation 
the asymptotics 

- At'V -1 

p hT (K, s) ~ — ——- —. — r-j (3.2) 

that we can rewrite in the form 

^ (K>8) ~ r(fe,r)|«h + ^ ^ r(/l ' r)= A^' (3 ' 3) 
Choosing r(/i, r) = 1 (it suffices to choose r(h,r) — >■ 1) we get 

^,r( K > s) ->■ p (k, s) = , (3.4) 

+ s p 

We honour by the name scaling relation our condition 

fih a 



At/ 3 



1. (3.5) 



Via r = (/i/A)/^") 1 ^ we can eliminate the parameter r, apply inverse 
Laplace transform to (3.2), fix k and send h — > 0. So, by the continuity 
theorem (for the Fourier transform of a probability distribution, see Feller 
(1971), . we can identify 



Po,o(^ s ) 



\K,\ a + S? 



as the Fourier-Laplace solution u(k, s) of the space-time fractional Cauchy 
problem (for x G H, t > 0) 

t D,f u(x, i) = s £> a u(x, i) , u(x, 0) = <J(x) , < a < 2 , < (3 < 1 . (3.6) 
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Here, for < ft < 1, we denote by t-Df the regularized fractional differential 
operator, see Gorenflo and Mainardi (1997), according to 

tD?g(t)= t Df l \g(t)-g(0)] (3.7) 

with the Riemann-Liouville fractional differential operator 

tD e g(t) ~)ra-f))dt h (t-t')f' ' (3.8) 

l s9 («), 0=1. 

Hence, in longscript: 

1 d /" 0(f) dr ^(o)r^ 



< 8 < 1 



^i 3 g(t) = { T fi- /3) rft y ° (t " tr r(1 " # ' (3-9) 

If </(t) exists we can write 

and the regularized fractional derivative coincides with the form introduced 
by Caputo, see Caputo and Mainardi (1971), Gorenflo and Mainardi (1997), 
Podlubny (1999), henceforth referred to as the Caputo derivative. Observe 
that in the special case ft — 1 the two fractional derivatives t D^ g(t) and 
tD 13 g(t) coincide, both then being equal to g'(t). 

The Riesz operator is a pseudo-differential operator according to 

j5ff = -\K\ a /(k) , welt, (3.10) 

compare Samko, Kilbas and Marichev (1993) and Rubin (1996). It has the 
Fourier symbol — 

In the transform domain (3.6) means 

s p - 1 u(K,s)-s f> - 1 = -H a M(/t, s) 

hence 

- - 1 

s) = — (3-11) 
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and looking back at (3,4 ) we see: u(k,s) = p 00 (k,s). Thus, under the 
scaling relation (3.5), the Fourier-Laplace solution of the CTRW integral 
equation (2.1) converges to the Fourier-Laplace solution of the space- 
time fractional Cauchy problem (3.6), and we conclude that the sojourn 
probability of the CTRW converges weakly (or "in law") to the solution 
of the Cauchy problem for the space-time fractional diffusion equation for 
every fixed t > 0. Later in this paper we will present another way of passing 
to the diffusion limit, a way in which by decoupling the transitions in time 
and in space we circumvent doubts on the correctness of the transition. 

For a comprehensive study of integral representations of the solution to 
the Cauchy problem (3.6) we recommend the paper by Mainardi, Luchko 
and Pagnini (2001). 

Subdiffusive and Superdiffusive Behaviour 

With regard to the parameters a and /3 in equation (3.6) we single out the 
cases (i), (ii) and (iii) by attributing names to them. 

(i) a — 1 , (5 = 1: normal or Gaussian diffusion (according to ^ = f^f), 

(ii) a — 2 , < j3 < 1: time-fractional diffusion, 

(iii) < a < 2 , (5 = 1: space-fractional diffusion. 

Let us now consider (compare Gorenflo nnd Abdel-Rehim (2004)) the 
equation (3.6) and its solution (3.11) in transform space that describe the 
evolution of the sojourn probability density u(x, t) of a wandering particle 
starting in the origin x = at the initial instant t — 0. We call this 
behaviour subdiffusive if the variance 



behaves for t — > oo like a power t 1 with < 7 < 1 , normal if 7 = 1, 
superdiffusive if 7 > 1 or if this variance is infinite for positive t. Using the 
fact that by Fourier transform theory 



and writing the right hand side of (3.11) as an infinite series in powers of 
k/s (convergent for s > 1) we find by termwise Laplace inversion 





u(k, t) = 1 



- + ••• = E p (-\ K \ a tP) 



r(i + /3) r(i + 2/3) 
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from which for t > we obtain, for all < (3 < 1 the result, 



2^/r(l + /3), if a = 2, 
oo , if < a < 2 . 



In the special case (i) of Gaussian diffusion this reduces to (<j(t)) 2 = 2t. 



4. Thinning (Rarefaction) of a Renewal 
Process under Power Law Regime 



We are going to give an account of the essentials of thinning a renewal 
process with power law waiting times, thereby leaning on the presentation 
by Gnedenko and Kovalenko (1968) but for reasons of transparency not 
decorating the power functions by slowly varying functions. Compare also 
Mainardi, Gorenflo and Scalas (2004) and Gorenflo and Mainardi (2008). 

Again (as in Section 2) with the t n in strictly increasing order, the time 
instants of a renewal process, = to < t\ < t2 < ■ ■ ■ , with iid waiting times 
Tk = tk ~ tk-i (generically denoted by T), thinning (or rarefaction) means 
that for each positive index k a decision is made: the event happening in the 
instant tk is deleted with probability p (where < p < 1) or is maintained 
with probability q — 1 — p. This procedure produces a thinned (or rarefied) 
renewal process, namely one with fewer events. Of particular interest for 
us is the case of q near zero which results in very few events in a moderate 
span of time. To compensate for this loss (wanting to keep a moderate 
number of events in a moderate span of time) we change the unit of time 
which amounts to multiply the (numerical value of) the waiting time with 
a positive factor r so that we get waiting times rTk and instants rtk in the 
rescaled process. Loosely speaking, it is our intention to dispose on r in 
relation to the rarefaction factor q in such way that for very small q in some 
sense the "average" number of events per unit of time remains unchanged. 
We will make these considerations precise in an asymptotic sense. 

Denoting by F(t) = P(T < t) the probability distribution function of 
the original waiting time T, by f(t) its density (generically this density is a 



generalized function represented by a measure) so that F(t) = / f(t') dt', 

Jo 

and analogously for the functions F k {t) and /&(£), the distribution and 
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density, respectively, of the sum of k waiting times, we have recursively 

fi(t) = f(t),f k (t)=[ t f k _ l (t)dF(t>) for k>2. (4.1) 

Jo 

Observing that after a maintained event of the original process the next one 
is kept with probability p but dropped with probability q in favour of the 
second-next with probability pq and, generally n — 1 events are dropped in 
favour of the n-th next with probability p n ~ 1 q, we get for the waiting time 
density of the thinned process the formula 

oo 

gM^Q^Mt). (4.2) 

n=l 

With the modified waiting time tT we have P(tT < t) — P(T < t/r) — 
F(t/r), hence the density f(t/r)/r, and analogously for the density of the 
sum of n waiting times f n (t/r)/r. The density of the waiting time of the 
sum of n waiting times of the rescaled (and thinned) process now turns out 

as 

oo 

g q ,r(t) = J2<lP n - 1 fn(tM/T. (4.3) 
n=l 

In the Laplace domain we have f n (s) = (/(s)) n , hence (using p—l — q) 

Us) = E ^ = : f iS \ 7( , ■ ( 4 - 4 ) 

n=l 1 -(!-?)/(«) 

By rescaling we get 

oo 

ft,r(s)=E?P B - 1 (7(r a )) 

n=l 

Being interested in stronger and stronger thinning (infinite thinning) let 
us consider a scale of processes with the parameters q of thinning and r of 
rescaling tending to zero under a scaling relation q = q(r) yet to be specified. 

Let us consider two cases for the (original) waiting time distribution, 
namely as in Lemma 2 of Section 3 case (A) of a finite mean waiting time 
and case (B) of a power law waiting time. We assume 

POO 

A:= / t' f(t')dt' <oo (A), setting = 1, (4.6A) 
Jo 



qf{rs)^ 

1 - (1 - q)J{rs) 



(4.5) 
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or 



POO 

(t) = / f(t') dt' ~ -t~ p for t ->■ oo with < (5 < 1 . (4.65) 
it P 



In case (B) we set 

A 



C7T 



+ 1) sin(/37r) ' 

By Lemma 2 of Section 3 we have f(s) = 1 — As' 3 + o(s^) for < s — > 0. 

Passing now to the limit q — > of infinite thinning under the scaling 
relation 

q = \t? (4.7) 

for fixed s the Laplace transform (4.5) of the rescaled density g qtT (t) of the 
thinned process tends to g(s) = 1/(1 + s 13 ) corresponding to the Mittag- 
Leffler density 

g(t) = -±E p (-lfi) = <ff L (t). (4.8) 

Thus, the thinned process converges weakly to the Mittag-Leffier renewal 
process described in Mainardi, Gorenflo and Scalas (2004) (called fractional 
Poisson process in Laskin (2003)) which in the special case (3 = 1 reduces 
to the Poisson process. In this sense the Mittag-Leffler renewal process is 
asymptotically universal for power law renewal processes. 



5. Mittag-Leffler Waiting Time and Space- 
Time Fractional Diffusion 

Let us sketch how, under the power law assumptions of Lemma 1 and 
Lemma 2, the simultaneous passage to the limit in inversion of the Fourier 
and Laplace transforms in (3.4) can be circumvented. Leaning on our 
presentations in Mainardi et al. (2000), Gorenflo et al. (2001) and Scalas, 
Gorenflo and Mainardi (2004), we introduce a memory function H(t) via 
which we will arrive at an evolutionary integral equation for the sojourn 
probability density p(x,t). For illustration we will soon consider a few 
special choices for this function. By rescaling and respeeding the process 
in time t and passing to an appropriate limit we will get a time-fractional 
evolution equation for p(x, t) (in fact a time-fractional generalization of 
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the Kolmogorov-Feller equation) that arises also by direct insertion of the 
Mittag-Lefner waiting time density into the CTRW integral equation (2.1) 
as we can see in Hilfer and Anton (1995). Via a second respeeding, obtained 
by rescaling the spatial variable x , we will arrive at the Cauchy problem 
(3.6) for space-time fractional diffusion. We keep the notations of Sections 
2 and 3. 

First, we introduce in the Laplace domain the auxiliary function 

S(f)(s) 0(s) 

and see by trivial calculation that (2.4) is equivalent to 

H(s) sp(n, s) - 1 = [w(k) - l]p(«, s) , (5.2) 

meaning in the space-time domain the generalized Kolmogorov-Feller 
equation 

rt Q r+oc 

H(t-t')—p(x,t')dt' = -p(x,t) + w(x - x')p(x',t)dx' , (5.3) 

JO & J-OD 

with p(x, 0) = S(x). 

Note that (5.1) can be inverted to 

?(*) = - — Wt- ( 5 - 4 ) 

1 + sH(s) 

We may play with equation (5.1), trying special choices for H(s) to obtain 
via (5.4) meaningful waiting time densities (j>(t). Or we chose H(t) hoping 
again to get via (5.4) a meaningful density <f)(t). In accordance with our 
inclination towards power laws let us take H(s) = s 13 ' 1 and distinguish the 
cases (i)/3 = 1 and (ii) < ft < 1. 
Case (i) yields 

H(s) = 1 , H(t) = S(t) , 0(s) = ^ , 0(t) = exp (-0 , 

namely the exponential waiting time density, and (5.2) reduces to 

sp(K, s) — 1 = [w(k) — 1] p(k, s) , (5.5) 



Rudolf GORENFLO 



15 



in the space-time domain the classical Kolmogorov- Feller equation 



H{t) = f -^y =: Hf L (t), 0( S ) = cf>(t) = ~E fi {-lP) = <Pff\t), 



namely the Mittag-Leffler waiting time density introduced in Section 1 by 
formula (1.4). With the Caputo derivative operator t D^ of (3.7) we get the 
Cauchy problem 



with p(x, 0) = 5(x). 

Remark: Because of (3.7) the equations (5.7) and (5.6) coincide in the 
special case (5 = 1. 

Rescaling and Respeeding 

Let us now manipulate the generalized Kolmogorov-Feller equation (5.3) by 
working on it in the Laplace domain via (5.2). 

Rescaling time means: With a positive scaling factor r (intended to be 
small) we replace the waiting time T by tT . This amounts to replacing 
the unit 1 of time by 1/r, and if r << 1 then in the rescaled process there 
will happen very many jumps in a moderate span of time (instead of the 
original moderate number in a moderate span of time). The rescaled waiting 
time density and its corresponding Laplace transform are (fi T (t) = (p(t/r)/r, 
(f) T (s) = <P(ts). Furthermore: 




(5.6) 



with p(x, 0) = d~(x) 
Case (ii) yields 




(5.7) 




1 -Ms) _ 1 -j>{rs) 
S(f) T (s) S(p(rs) 



, hence T (s) 



l + sH T (s) ' 



1 



(5.8) 



and (5.2) goes over into 



H T (s) sp T (n, s) - 1 = [w(k) - 1] p t (k, s) . 



(5.9) 
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Remark: Note that in this Section 5 the position and meaning of the 
indices attached to the generic density p are convenient but different from 
those in Section 3. 

Respeeding the process means multiplying the left hand side (actually 

d 

—p(x,t') of equation (5.3) by a positive factor 1/a, or equivalently its 

right hand side by a positive factor a. We honour the number a by the 
name respeeding factor, a > 1 means acceleration, a < 1 deceleration. In 
the Fourier-Laplace domain the rescaled and respeeded CTRW process then 
assumes the form, analogous to (5.2) and (5.9), 

HtM sp TA (n, s) - 1 = a [w(k) - 1] P T)0 (k, s) , (5.10) 
with _ _ 

a as<f)(rs) 

What is the effect of such combined rescaling and respeeding"! We find 

4>ra{s) = 1 = a * (T3 l , (5.11) 

l + sH T>a (s) 1 - (1 - a) 0(r S ) 

and are now in the position to address the Asymptotic universality of 
the Mittag-Leffler waiting time density. 

Using Lemma 2 with rs in place of s and taking 

a = \r fS , (5.12) 

fixing s as required by the continuity theorem of probability for Laplace 
transforms, the asymptotics 4>{rs) = 1 — A(rs) /3 ) + o((rs) /3 ) for r — > 
implies 

% M Ar^l-A(r S )Q + o((r S n] 1 7ml r , 

corresponding to the Mittag-Leffler density 

<t>f L {t) = -jE,{-t?) . 
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Observe that the parameter A does not appear in the limit 1/(1 + s 13 ). 
We can make it reappear by choosing the respeeding factor r 13 in place of 
Xt 13 . In fact: 

Formula (5.13) says that the general density <f)(t) with power law 
asymptotics as in Lemma 2 is gradually deformed into the Mittag-Lefhcr 
waiting time density <^ L {t). It means that with larger and larger unit 
of time (by sending r — > 0) and stronger and stronger deceleration (by 
a = Xt 13 ) as described our process becomes indistinguishable from one with 
Mittag-Leffler waiting time (the probability distribution of jumps remaining 
unchanged). Likewise a pure renewal process with asymptotic power law 
density becomes indistinguishable from the one with Mittag-Leffler waiting 
time (the fractional generalization of the Poisson process by Laskin (2003) 
and Mainardi, Gorenflo and Scalas (2004)). In fact, we can consider the 
pure renewal process as a CTRW with jump density w(x) = 5(x — 1), the 
position of the wandering particle representing the counting number (the 
number of events up to and including the instant t). 

Remark: It is instructive to look at the effect of combined rescaling 
and respeeding on the Mittag-Leffler density 4>f L {t) itself which by (1.6) 

also obeys the asymptotic conditions of Lemma 2. We have 4>^ L (s) = 
1 - s p + o(s p ) for s -+ 0, and with (5.11) we find the relation 

Uf L \ ( s ) = 0f L (Ts/a 1/l3 ) for all r > , a > , (5.14) 

expressing the self-similarity of the Mittag-Leffler density. In particular we 
have the formulas of self-similarity and invariance 

« £ ) T>T , (f) = < L (0 for all r>0, (5.15) 

telling us that the Mittag-Leffler density is invariant under the transforma- 
tion (5.11) with the respeeding factor a = t 13 in place of a = Xt 13 . 

Diffusion Limit in Space 

In addition to rescaling time we now rescale also the spatial variable x, 
by replacing the jumps X by jumps hX with positive scaling factor h, 
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intended to be small. The rescaled jump density turns out as Wh{x) = 
w(x/h)/h, corresponding to Whi^) = w(kK). Starting from the Fourier- 
Laplace representation (5.2) of our CTRW with general waiting time density, 
we accelerate the spatially rescaled process by the respeeding factor l/(fih a ) 
with ji > and arrive (using as new dependent variable) at the equation 



H(s) 



s) - 1 



w{kK) 



fih a 



(5.16) 



Assuming now the power law assumptions of Lemma 1 satisfied, fixing k 
as required by the continuity theorem for the Fourier transform, we get 
(w(nh) — l)/(/x/i a ) — > —\n\ a for h — > 0, and writing u in place of q , in the 
limit 

H(s) su(k,s)-1 = -\n\ a u(K, s) . (5.17) 
By Fourier inversion we get 

H{s) [su(x, s) - S(x)] = x Dq u(x, s) , 

and then by Laplace inversion in the space-time domain the Cauchy problem 



1 d 

H(t - t') j-u{x, t) = X D* u(x, t) , u(x, 0) = S(x) , < a < 2 . (5.18) 



As in Section 3 we mean by x Dq the Riesz pseudo-differential operator with 
Fourier symbol — \k\ according to formula (3.10). 

Finally, inserting into (5.18) the Mittag-Leffler memory function 



6(t) ;i£P = l, 



;if <0<1, 



we recover the Cauchy problem (3.6) for the space-time fractional diffusion 
equation, namely 



t D^u(x,t) = x D£u(x,t) , u(x,0) = 5(x) ,0<«<2,0</3<l .(5.19) 



Comments: In this Section 5 we have split the passage to the limit into 
a temporal one r — > and a spatial one h — > 0. In Section 3 by the 
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well-scaled (combined) passage to the limit we have avoided the concept 
of respeeding but have had to eliminate the parameter r via the scaling 
relation (3.5). But where has the scaling relation (3.5) gone? We rediscover 
it as hidden in the deceleration (5.10) with (5.12) and the compensating 
acceleration (5.16). We recommend to compare this splitting with the more 
abstract technique of triangular arrays applied by Meerschaert and Schemer 
(2004) and (2008). We think that our method offers intuitive insight into 
the meaning of passing to the diffusion limit. 

6. Time-Fractional Drift and Subordination 

In Section 2 we have seen in which way a CTRW is subordinated to 
a renewal process , and in Section 5 we have worked out the effect of 
subordination under the Mittag-Leffler renewal process, see equation (5.7). 
For the following considerations we hint to the references Mainardi, Gorenflo 
and Scalas (2004), Meerschaert et al (2002), Mainardi, Pagnini and Gorenflo 
(2003), Mainardi, Luchko and Pagnini (2001), Gorenflo, Mainardi and 
Vivoli (2007), Gorenflo and Mainardi (2008), and again to the papers by 
Meerschaert and Scheffler (2004) and (2008). 

Our aim is to pass to a meaningful limit of the Mittag-Leffler renewal 
process whose waiting density is 4>f L {t) under the restriction < /3 < 1, 
meaning exclusion of the limiting case j3 — 1 of exponential waiting time. 

Keeping our earlier notations, in particular Tk for the waiting times, tk 
for the jump instants, and denoting by N = N(t) = max {k\tk < t} the 
number of renewal or jump events up to instant t (the "counting number") 
we have for a general renewal process the probability 



Finding it convenient to embed the formalism of renewal into the CTRW 
formalism we introduce a pseudo-spatial variable r taking the values of 
the counting number N and denote the sojourn probability by q(r,t) (with 
r > 0, t > 0), by v(r) the jump density. Taking (because N runs through 
the non-negative integers) for the jump- width of the "random walk" so 
generated the constant 1 we have v(r) = S(r — 1) and v(k) = exp (in), and 
(5.2) yields 



P(N(t) = k) = P(t k <t,t k+1 >t). 



(6.1) 



H(s) sq(n, s) — 1 = \v(k) — 1] q(n, s) . 



(6.2) 
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With the waiting time density (fi^ L (t) and correspondingly 

1 



lML, 



s = 



1 + 



H(s) 



(6.2) goes over in 



J3-1 



sq(n,, s) — 1 = [exp {in) — 1] q(n, s) . 



(6.3) 



As we have done in Section 5 for the general CTRW we now ask what 
happens when we pass to the diffusion limit in "space" for the Mittag- 
Leffler renewal process in the CTRW formalism. Multiplying the jumps 
by a positive scaling factor 5, decorating q by such index, replacing v(k) 
by v(k,5) according to vs{r) = v(r/5)/5, fixing k,, finally accelerating by 
applying the factor 5~ l to the right hand side we obtain the equation 



sq s {K, s) 



5 1 [exp (inS) — 1] q s (n, s) . 



and 5 — > yields the equation 

= iK%{K,s) (6.4) 
which implies 

f («, a ) = -i- 1 -. (6.5) 

We remark that an analogous limit can by proper scaling also be obtained 
directly from the generic power law renewal process that we have discussed 
in Section 4 on thinning. 

Note that (6.4) corresponds to the Cauchy problem for the positively 
oriented fractional drift equation with the Caputo derivative operator t D@ , 

t D%(r,t) = -^q (r,t) . (6.6) 

Without inverting the transforms in (6.5) we can recognize the self- similarity 
of the function q (r,t). With any positive constants a and b, generic 
functions / and g we have the correspondence of f(ax) to /(/c/a)/a 



J3-1 



sq (K, s) 



Rudolf GORENFLO 



21 



and the correspondence of g(bt) to g(s/b)/b, hence via (6.5) qo{ax,bt) = 
b" 13 q (ax/bP ,t), and with Qq[x) = qo(x, 1) and the similarity variable x/t 13 
we obtain 

q (x,t)=t-PQ (x/t?). (6.7) 
Fourier inversion of (6.5 ) gives 

( s 13 ' 1 exp (-rs 13 ) , forr > 
, for r < 



Using the fact that exp (— s 13 ) is the Laplace transform of the extreme 
positive-oriented unilateral stable density of order (3, namely of Lp(t), we 

^ qo(r,t) = r-W J 1 -* L? {tr-W) , (6.9) 

with the Riemann-Liouville fractional integration 

t P git) = ^- J\t - ty- 1 g{t') dt', 7 > . 

This solution can be expressed in alternative ways for which we refer to the 
references cited at begin of this Section. We have, for t > 0, 

qo(r,t) = y-i-W L f(tr-W)=t-eM fi ( r t-f l ), (6.10) 

with the M- Wright function 

oo ( \n i 00 ( \n—l 

for whose properties and use we refer to Mainardi, Mura and Pagnini (2010). 

We recognize the function qo{r,t) as the sojourn probability density of 
the directing process (the subordinator) for space-time fractional diffusion 
(3.6), according to Meerschaert et al. (2002), Mainardi, Pagnini and 
Gorenflo (2003), Gorenfio, Mainardi and Vivoli (2007). 

Let us use the notations of Gorenflo, Mainardi and Vivoli (2007). There 
it is outlined how the random position x(t) of the wandering particle can 
be expressed via an operational time t* (the properly scaled limit of the 
counting number of the renewal process to which he CTRW is subordinated) 
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from which the physical time t = t(t*) is generated by a positive-oriented 
stable extremal process of order /3 whereas the spatial position is generated 
by a stochastic process y = y(t*), stable of order a. Then we have the 
process x = y(t*(t)) in parametric representation 



In practice of simulation this representation can be used to produce "flash- 
light shots" of a set of points in the (t, x) plane showing a sequence 
of positions of a wandering particle. One only needs routines for 
generating random numbers from the relevant stable probability laws. As 
a consequence of infinite divisibility and Markovianity of these laws the 
sequence of points so produced constitutes a sequence of true particle 
positions (possible positions of a really happening process), see Gorenflo, 
Mainardi and Vivoli (2007). 

However, it is more common, see Feller (1971), to treat subordination 
directly in the form x = x(t) = y(t*(t)), namely to produce first the process 
t* = t*(i) of generating the operational time t* from the physical time t. 

The processes t = t*(t) and t* = t*(t) are inverse to each other, but 
if < j3 < 1 the latter is neither Markovian nor infinitely divisible. Its 
probability density function is go(^*, t), evolving in physical time t and given 
by formula (6.10). Sometimes the process is called a Mittag-Leffler process 
(see Meerschaert et al. (2002) and search for further details also in several 
chapters of Feller's famous book). The motivation for calling the process 
after Mittag-Leffler lies in its various relation to the Mittag-Leffler function. 
The inverse Laplace transform of (6.5) is q (K,,t) = Ep(mt) (a Mittag- 
Leffler function with imaginary argument) and the M- Wright Mp has by 
Fourier and Laplace transform two more connections to the Mittag-Leffler 
function, namely, see Mainardi, Mura and Pagnini (2010). For < /3 < 1 
we have Mp(s) = Ep(—s) and the Fourier correspondence of Mg(|a;|) to 
2E 2( g(— k 2 ). Furthermore, see Bondesson, Kristiansen and Steutel (1996) 
and Meerschaert and Scheffler (2004), the Laplace- type integral 



* = t(U) 



x = y(U) . 




From (6.5) follows by Fourier inversion 




(6.11) 
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In the case < /3 < 1 the probabilty law of the process t* = £*(£) is for no 
positive t infinitely divisible, see Bondesson, Kristiansen and Steutel (1996) 
and Meerschaert and Scheffler (2004). This process must not be confused 
with the process that is called Mittag-Leffler process by Pillai (1990) and 
that is obtained via the infinite divisibility of the Mittag-Leffler distribution 
whose density is <pf L = --^E^-t 13 ). 

We display the resulting subordination formula (compare Meerschaert 
et al. (2002) and Gorenflo, Mainardi and Vivoli (2007)) for the solution 
u(x,t) of the Cauchy problem (3.6) (in formula (5.19) repeated): 



where the symmetric stable density f a (x,r) has the Fourier transform 
/ a («, t) exp (— r|/t| a ) and the variable r represents the operational time t*. 

Let us finally observe that in the excluded limiting case /3 — 1 by again 
identifying r with t*, equation (6.8) leads to the result qo{t*,t) = S(t — i*), 
hence t — t* which means that in the case of exponential waiting time the 
physical time and the operational time coincide. 



We have discussed some (essentially two) ways of passing to the diffusion 
limit from continuous time random walk with power laws in time and 
in space (for transparency of presentation spatially symmetric and one- 
dimensional), namely 

(i) what we call well-scaled passage to the limit where the rescalings of time 
and space are carried out in a combined way, 

(ii) carrying out the passages separately in time and in space. 

The limit in time, by aid of a convenient memory function, leads to the 
Mittag-Leffler waiting time renewal process or fractional Poisson process, 
and the Mittag-Leffler function becomes essential for description of long- 
time behaviour (of renewal processes and of CTRW) whereas for the wide- 
space view stable distributions take the role. In the time domain there are 
two passages to the limit. The first one leads to an extremal stable density 
evolving in time, the other one by condensing the corresponding counting 
process to smaller and smaller counting-step-size leading to a Mittag-Leffler 
process as the subordinator of the continuous time random walk. 




(6.12) 



7. Conclusions 
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Actually we have obtained this subordinates: by another splitting: by 
passing first to the Fractional Poisson process and from this then to the 
subordinator. But this additional splitting can be avoided. For our way of 
analyzing the transition to the limit in power law renewal processes we have 
got inspiration from studying the theory of thinning such processes, and we 
have discovered the important analogy of a limit formula in the Laplace 
domain. 

Quite generally, for performing the necessary investigations in CTRW 
theory the Fourier-Laplace domain is the most convenient operational 
playground since a long time. To make well visible the basic ideas we have 
avoided measure-theoretic and functional- analytic terminology, hoping so 
to be not too difficult for people not so well trained in these fields. 
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